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We provide a parton construction of wavefunctions and effective field theories for fractional Chern insulators. 
We also analyze a strong coupling expansion in lattice gauge theory that enables us to reliably map the parton 
gauge theory onto the microsopic Hamiltonian. We show that this strong coupling expansion is useful because 
of a special hierarchy of energy scales in fractional quantum Hall physics. Our procedure is illustrated using the 
Hofstadter model and then applied to bosons at 1/2 filling and fermions at 1/3 filling in a checkerboard lattice 
model recently studied numerically. Because our construction provides a more or less unique mapping from 
microscopic model to effective parton description, we obtain wavefunctions in the same phase as the observed 
fractional Chern insulators without tuning any continuous parameters. 



INTRODUCTION 

The discovery of the fractional quantum Hall effect in lay- 
ered semiconductor devices in high magnetic field [1] opened 
our eyes to a rich world of truly quantum mechanical phases of 
matter that exist at zero temperature. Quantum entanglement 
plays a crucial role in defining and understanding topological 
phases of matter like the fractional quantum Hall liquids. In- 
stead of long range order and symmetry breaking, we should 
study the pattern of long range entanglement in such phases 
[2, 3]). One manifestation of the necessity of quantum entan- 
glement is the inability of product or mean-field like wave- 
functions to capture, even qualitatively, the physics of such 
a phase. Their potential ability to function as quantum com- 
puters [4], as well as their natural robustness to local deco- 
herence, are also consequences of the presence of long range 
entanglement in the fractional quantum Hall fluids. 

Despite the beauty and ruggedness of these topological liq- 
uids, there are still many practical challenges as we attempt to 
observe non-Abelian particles in nature and construct the first 
scalable quantum computer. The zero-temperature robustness 
of these topological fluids is eventually lost at finite temper- 
ature, and since all experiments are carried out at finite tem- 
perature, we are forced into practical questions. Exactly how 
big is the gap to excitations? Just how easy is it to implement 
the non-local operations that detect non-Abelian excitations or 
perform useful quantum computations? Moreover, given the 
conceptual and practical importance of topological phases of 
matter, it is very interesting to ask where else in nature such 
topological phases may be lurking. 

Motivated by these questions and others, there has been an 
explosion of interest in new models that mimic the physics of 
the quantum Hall effect. In the Hall effect, the basic starting 
point is Landau levels, perfectly flat bands that exist in a uni- 
form magnetic field. When Landau levels are totally filled we 
find the integer quantum Hall effect, and when Landau levels 
are partially filled, there exists a large degeneracy in the non- 
interacting limit. The inclusion of interactions resolves the 
degeneracy and produces an incompressible topological fluid 
[5], a fractional quantum Hall fluid. The physics of Landau 
levels may be conjured anew on the lattice by studying tight 



binding models that break time reversal symmetry [6, 7]. Such 
models can be tuned to give very flat bands which simultane- 
ously possess a non-trivial Chern number. Filled bands with 
non-trivial Chern number are in the same universality class as 
integer quantum Hall states, and the phenomenology is largely 
transferable. Thus it is reasonable to guess that if we were able 
to partially fill a flat Chern band and expose the electrons to 
interactions of a magnitude much greater than the bandwidth, 
then the physics of the fractional quantum Hall effect should 
also be visible (see Ref. [8] for a discussion of this point). This 
beautiful idea, which is a striking example of universality in 
physics, has recently been numerically confirmed in a flurry of 
activity [9-13]. (Previous work, motivated by the possibility 
of cold atoms realizations of FQHE includes [14, 15].) 

From the perspective of the questions raised above, these 
fractionalized Chern insulators are interesting in that they pro- 
vide a new realization of quantum Hall physics with the po- 
tential for much higher energy gaps and the possibility of new 
methods for the detection anyons. Particularly interesting is 
the idea of realizing new high temperature non-Abelian topo- 
logical liquids, on which we comment later. As an important 
early step towards the observation and manipulation of non- 
Abelian anyons in fractional Chern insulators, we must es- 
tablish a basic theoretical framework in which to understand 
numerical and experimental results. We have in mind the pow- 
erful formalism in fractional quantum Hall physics known as 
the slave-particle or parton approach that provides wavefunc- 
tions and low energy effective theories for fractional Hall liq- 
uids (e.g. [3, 16-18]). In practice, the parton approach to this 
problem represents a rationalization for and a massive gen- 
eralization of Laughlin's original wavefunction [19] for the 
f = 1/3 plateau. We are interested in transferring this tech- 
nology to the new setting of fractional Chern insulators. An 
interesting recent attempt to make this transition and provide 
wavefunctions for fractional Chern insulators can be found in 
Ref [20]. Parton wavefunctions similar to those we provide 
have been used in variational studies in Ref. [21]. Ref [22] 
also uses technology from quantum Hall by embedding the 
Chern band problem into a larger Landau level. Ref. [23] ex- 
tends a Woo structure of quantum Hall states to Chern bands. 

In this paper we achieve this aim by formulating a parton 
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description of the recently observed Abelian fractional Chern 
insulators. This parton description provides model wavefunc- 
tions and low energy effective theories that demonstrate the 
universality of fractional quantum Hall physics. The core of 
the parton formalism is a mapping of the original electronic 
system to an effective description in terms of fractionally 
charged "partons" coupled to emergent gauge fields. Analysis 
of a parton model always requires analysis of a gauge theory, 
and in the case of fractional Chern insulators, we deal with 
lattice gauge theory. In addition to the basic story involving 
wavefunctions, we formulate and carry out a strong coupling 
expansion of certain lattice gauge theories. This strong cou- 
pling analysis' yields two immediate benefits: new insights 
into the dynamics of some strongly coupled gauge theories, 
and new clues leading to fractional Chern insulators, includ- 
ing novel non- Abelian states. In fact, because the strong cou- 
pling expansion provides an almost unique mapping from the 
microscopic Hamiltonian to the lattice gauge theory, our con- 
struction can be regarded as a zero parameter model of the 
original microscopic physics in terms of the gauge degrees of 
freedom. Because the gauge theory dynamics can be reliably 
analyzed (thanks to the Chern-Simons term), our construction 
provides strong evidence that the system enters a fractional 
Chern insulating phase. 

As we explain in more detail below, the main conceptual 
obstacle to the construction of a parton theory for fractional 
Chern insulators is the question of the parton band structure. 
In the case of Landau levels, the degeneracy of each Landau 
level naturally depends on the charge of the particle moving 
in the magnetic field. Consider the case of = 1 /3. Breaking 
the electron into three "colors" of charge 1/3 partons leads to 
a partonic Landau level with reduced degeneracy, and since 
the number of partons of each color is the same as the number 
of electrons, the partons are immediately able to completely 
fill their Landau level. This picture leads to an easy accounting 
of the Hall conductivity, as we have three colors of partons 
each contributing a Hall conductivity of ''^^-jp— for a total Hall 
conductivity of | ^ . Vaezi [26] has noticed that the story for 
fractional Chern insulators is not as simple, since it is less 
straightforward to reduce the size of the partonic bands, but 
the state Vaezi constructs is not obviously gauge invariant and 
appears to give zero after projection. We address this issue in 
our construction by expanding the size of the unit cell seen by 
the partons in a way that remains invisible to the microscopic 
electrons. This discrete set of choices is the only freedom in 
the parton band structure that is needed to open a gap in the 
single particle spectrum. 

The remainder of the paper is organized as follows. First, 
we describe our construction in the context of a simple lat- 
tice model of the fractional quantum Hall effect, namely the 



Deriving the dynamics of liadronic objects using a strong coupling expan- 
sion in lattice gauge tlieory dates back at least to [24] and is used in the 
same spiiit in e.g. [25]. 



partially-filled Hofstadter model. In the same section we dis- 
cuss a strong coupling expansion for the parton gauge theory 
that, owing to the special physics of the fractional quantum 
Hall effect, leads to an unusually complete physical justifi- 
cation for the parton approach. Second, we apply our con- 
struction to a case that makes contact with recent numerical 
calculations, fractionally filled fermions and fractionally filled 
bosons moving in a checkerboard model. Finally, we summa- 
rize our results and indicate some exciting directions we are 
currently exploring. 

Note: while this paper was in preparation, we learned of a 
paper by Lu and Ran [27] which also develops a parton con- 
struction for fractional Chern insulators. The overlap between 
their work and ours is significant; our discussion of the strong 
coupling expansion is a substantial difference. Their work is 
a general discussion of a variety of parton wavefunctions in- 
cluding time reversal symmetric insulators, while we focus 
on the story of the strong coupling expansion and the case of 
fractional Chern insulators. 



HOFSTADTER MODEL AND STRONG COUPLING 
EXPANSION 

We introduce our construction in the simple context of the 
Hofstadter model of particles hopping in a uniform magnetic 
field on a lattice [28]. The Hofstadter model is defined by 
taking a square lattice and placing 2tt/N flux through each 
plaquette of the lattice. For a fixed N, the Hofstadter model 
has a unit cell consisting of N sites and thus a band structure 
consisting of N bands. As TV — > oo, the lower bands of the 
Hofstadter model reproduce the continuum physics of Lan- 
dau levels. In terms of universal physics, the lowest band of 
the Hofstadter model has Chern number one and hence repro- 
duces all the low energy physics of the v = 1 integer quantum 
Hall state even at finite N. Similarly, a fractionally filled low- 
est Hofstadter band in the presence of interactions is a natural 
lattice regularizaton of a fractional quantum Hall fluid. We 
now consider the case of electrons moving in the lowest Hof- 
stadter band at filling 1 /3. 

Consider the amplitude for an electron to hop around one 
plaquette of the square lattice. This amplitude has a phase 
given by e^^*/^ which we interpret as an Aharanov-Bohm 
phase for an electron of charge e moving in a field of mag- 
nitude B = j^^^2 where a is the lattice spacing. Although 
the models discussed in [9-13] generally speak of quantum 
Hall type systems without magnetic fields, we see that we can 
always reinterpret phases among the hopping parameters in 
terms of strong lattice scale magnetic fields. This interpre- 
tation is useful because it suggests a way to obtain a sensi- 
ble parton band structure. If the electron fractionalizes into 
charge e/3 partons then these partons should naturally expe- 
rience a phase of e^'^'/'^'^^) when moving around a plaquette 
because of their reduced charge. Thus at the mean field level, 
before including the gauge fluctuations that glue the partons 
back together into electrons, the partons move in a modified 
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Hofstadter bandstructure with N replaced by 3 A^. This means 
the partons move in an enlarged unit cell relative to the elec- 
trons, and with one parton of each color per electron, each 
parton color can fill an entire band, leading to a gapped mean 
field ground state. 

Despite the apparent enlargement of the unit cell, the phys- 
ical electron does not experience an enlarged unit cell since 
it encircles 27r flux after moving around only N plaquettes. 
Equivalently, motion of the electron requires coordinated mo- 
tion of the partons, and counting the phase of 27r/3 per color 
accumulated by a parton around plaquettes, we find a total 
phase 2tt after summing over the three colors. The analogy to 
the usual continuum story [3, 17, 18] is, we hope, now clear. 
In particular, the electron wavefunction is obtained by taking 
the mean field wavefunction of the partons, three copies of 
the lowest 3 Hofstadter band, and projecting onto the color- 
neutral state. As we take iV — > cx) this procedure recovers 
Laughlin's wavefunction for = 1/3. 
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FIG. 1: A rendition of the checkerboard model showing the NNNN 
hoppings interpreted as "bridges" out of the plane connecting more 
distant sites. The red and blue dots are the two sublattices. Our 
lattice gauge theory version of this model associates an SU{3) link 
variable with every elementary hopping path shown here. Included in 
the sum over loops are loops running along these out of plane bridges 
representing the NNNN hoppings. 

The above procedure may seem somewhat ad hoc, but there 
is a more systematic way to relate the parton band structure to 
the electron band structures from which it emerges. First, we 
formalize the parton construction by writing the electron op- 
erator Cr at site r in terms of three partons fra carrying charge 

e/3 as Cr = frlfr2fr3 = ^ £0/^7 /ra /r/3 /r7 aS in Rcf. [17]. 

This expression is manifestly symmetric under SU{3) trans- 
formations of the fra, and since we can make such transfor- 
mations at each site r, we have an SU{3) gauge structure. 
The group SU{3) is the so-called high energy gauge group. 
The terminology is necessary because the mean field parton 
Hamiltonian may break some of the SU{3) symmetry, so the 
true low energy gauge group can be distinct from SU{3). 
Next, we write a Hamiltonian for the parton theory includ- 
ing the SU{3) gauge fluctuations. For every term in the par- 
ton hopping Hamiltonian we associate an SU{3) link variable 
V. This particular choice is very convenient for performing 



a strong coupling expansion, and can be visualized as allow- 
ing the partons to hop longer distances through lattice scale 
"bridges" or "wormholes" as shown in Fig. 1. Our notation is 
as follows: the lattice is defined by a graph (V,£), fermions 
fr live on the sites r G V and gauge fields Vrr' live on links 
rr' e £. There is a link in £ for every elementary hopping 
term in the parton Hamiltonian. The set of elementary loops 
on the lattice we denote C The (chromo)electric field conju- 
gate to Vrr' is E:^r' whcrc A is an adjoint index for SU{3). 
With these conventions, our Hamiltoruan is 

^ = - E Cf^Vrr'fr + h.C. + ^ E^, 
rr'^£ rr' ^£ 

-KY,^rl\{VrrA+h.c. (1) 

leC \rr'el / 

The first term is a parton hopping term, the second electric 
term favors small electric fields, and the third magnetic term 
favors smooth gauge configurations. E'^^, is the quadratic 
Casimir of SU (3) on the link rr' , so the eigenstates of the 
electric field term are labelled by irreducible representions of 
SU{?>). In this basis, the magnetic field term functions like a 
raising and lowering operator: it adds an electric field line in 
the fundamental along the path £. 

Note that the physical Hilbert space of the gauge theory 
may be formally larger than the electron Hilbert space due to 
the presence of arbitrarily highly excited gauge fields states 
or "glueballs". Thus the gauge theory technically describes 
electrons coupled to some high energy tower of bosonic states, 
however, these bosonic degrees of freedom are irrelevant for 
the low energy physics. We thus proceed to study the gauge 
system at strong coupling as a model of the low energy physics 
of electrons. 

Let us analyze this lattice gauge theory in a strong coupling 
expansion h ^ K, tl^, . We will need only the lowest order 
terms for our purposes here, but the expansion can be carried 
out systematically to higher orders where it is important to re- 
move disconnected terms [29]. The excitations that survive in 
the /i — > oo limit are the colorless fermionic "baryons" gen- 
erated by the operators Cr — frifr2fr3, that is, the original 
electrons. Away from h — oo, these electrons can hop via 
virtual fluctuations of their constituent partons. Consider a 
particular parton hopping amplitude and suppose we begin 
with an electron on site r. The parton hopping Hamiltonian 
moves one parton from site r to site r' with amplitude t-^ , but 
it also creates an electric field line connecting r to r'. Thus 
this intermediate state costs an energy of roughly h. We now 
proceed to hop a second parton with the resulting intermedi- 
ate state still containing electric fields. Finally, we hop the 
third parton, remove all electric fields, and return to a color- 
less state at r' . This process is illustrated in panel {A) of Fig. 
2. In effect, the electron moved from r to r' with an ampli- 
tude given by ^ {t^Y /K^. This argument applies to all 
the parton hopping amplitudes, so we have for every parton 
hopping amplitude a corresponding electronic hopping ampli- 
tude given by [t^Y /h? (with the conventional overall minus 
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sign). Furthermore, in our scheme the prefactor is the same in 
all cases (only one link ever has an electric field excited). Of 
course, there are other contributions to the electron amplitude 
from more complicated processes, but these are further sup- 
pressed by powers of /h. Note that in this analysis we focus 
on the interactions between two baryons and ignore possible 
many-body effects from the background density of baryons; 
this should be justified since such effects will require more 
complicated parton exchange patterns and will be further sup- 
pressed. 

>- time 
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FIG. 2: The two most important processes in the strong coupling 
expansion. Process (A) generates an electron kinetic term of order 
{t^)"^ /h^ . Process (B) generates repulsive electron interactions of 
Qxdsr{t^fK/h''. 

If we express all electron hoppings in terms of a particular 
matrix element ig, say a nearest neighbor hopping, then we 

have the remarkable result that -jif- — . In other 

*o \4 J 

words, the electron hoppings measured relative to a refer- 
ence hopping are literally the cubes of the parton hoppings 
measured relative to the equivalent parton reference hopping. 
Hence the parton bandstructure represents, in a very precise 
sense, the "cube root" of the electron band structure. It is im- 
portant that the mapping from electrons to partons is not quite 
unique, we always have discrete phase choice not set by the 
electron band structure. In the case of the Hofstadter model, 
this prescription is very simple since there is only one kind of 
hopping term, and so, up to the overall energy scale, we re- 
produce our simple guess that the parton Hamiltonian should 
have a third of the phase of the electron Hamiltonian through 
every plaquette. 

So far we have only mentioned electron kinetic terms. Elec- 
tron interactions are also generated. A simple process that 
generates such interactions is parton exchange between two 
neighboring electrons as shown in panel (B) of Fig. 2. This 
process involves not only parton hopping but also the mag- 
netic term to remove the final closed field line, and so it 



has an amplitude ^ K{t^)^/h'^ and is repulsive (an ex- 
tra minus sign comes the fermion exchange). In the strong 
coupling limit electron interactions are parametrically smaller 
than electron kinetic terms. Interactions are also short ranged 
since parton exchange over a distance x is supressed by fac- 
tors of (t^ /hY . Note also that the magnetic term in the gauge 
Hamiltonian generates ring-exchange terms for the electrons 
[25], where electrons hop around loops in a coordinated fash- 
ion. These terms are of order Kltf /hf and can be para- 
metrically smaller than and V^. Now the usual story at 
this point would be simple: the gauge field confines and we 
can integrate it out to produce a weakly coupled electronic 
Hamiltonian which can be analyzed in the usual way. We em- 
phasize that confinement is a very common situation in gauge 
theory, and deconfinement typically requires special circum- 
stances: high density, many flavors, broken time reversal, etc. 
Since strongly correlated electron systems typically need elec- 
tron interactions of the same order as electron kinetic terms, 
such a phase corresponds to the medium coupling region in 
the gauge theory h ^ ^ K. The usual hope is that once 
we make K and large enough, the gauge theory will flow 
in the infrared to a deconfined phase, but we emphasize that 
because the strong coupling expansion breaks down at inter- 
mediate coupling, we cannot use it to reliably predict elec- 
tronic models that realize such deconfined phases. Of course, 
it remains invaluable as a source of intuition. 

However, the case of fractional Chern insulators is special. 
This is because although the bare electron hoppings are much 
larger at large h than the electron interactions, the bandwidth 
of the nearly flat electronic band can be tuned to be quite 
small. We use the Hofstadter model as an example. Let the 
electronic gap be A"^ ^ f^, let the bandwidth of interest be 
w'^, and suppose it is tuned to be some small fraction F of the 
electron hopping (in the Hofstadter model, this can be accom- 
plished by making N large; more generally it is the goal of 
much bandstructure engineering [30, 31]). Then we want to 
implement the hierarchy of scales A"^ 3> C/"^ 3> . This can 
in fact be achieved within a region of parameter space where 
the strong coupling expansion is reliable, so long as F <C 1. 

All we need is /i ^ Vtf K and h <C \J~^^^- Of course, unless 
the electronic band is perfectly flat, we will eventually violate 
the hierarchy if we take h to infinity or K to zero while fixing 
tf. 

Based on the model systems studied in Refs. [9-13], a rea- 
sonable target Hamiltonian for realizing a fractional Chem in- 
sulator should have the following features: 

1. Large electronic band gap: A^ w"^, V^, 

2. Strong electron interactions: f/^ ^ u;"^, 

3. Small ring exchange terms. 

Regarding criterion 3, it is usually believed that ring exchange 
terms are helpful in the fight against confinement and we have 
no reason to suspect otherwise here. Nevertheless, the most 
conservative path is to limit ourselves as much as possible to 
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the minimal ingredients present in numerical studies that re- 
alize fractional Chern insulators.^ Of course, ring exchange 
terms will be generated under renormalization regardless of 
whether we include them in the microscopic Hamiltonian. We 
can realize all three criteria within the strong coupling expan- 
sion by establishing the hierarchy \J {tf K)/F ^ h ^ K > 
. This guarantees the electronic hierarchy A'^ W ^ w'^ 
plus small ring exchange terms. 

We are now in a privileged position. We have a precise and 
controlled map between the gauge theory and the microscopic 
electron theory via the strong coupling expansion. We also 
know the fate of the electron theory: it forms an incompress- 
ible phase with the same universal physics as the = 1/3 
fractional Hall fluid. Thus we also know the fate of the gauge 
theory: it must deconfine. To be fair, it is possible that the par- 
ticular interactions generated somehow always favor a crys- 
talline state, but we regard this as very unlikely. Thus decon- 
finement must occur even if tf and K are relatively small. As 
we said above, h cannot be infinite and K cannot be zero, 
but our partonic "molecules" can be relatively tightly bound 
in isolation and still deconfine when put together into a liq- 
uid. We note that there exist duality arguments that suggest 
that a lattice Chern-Simons term pushes the confinement tran- 
sition to /i = CXI [32]. A hint that this is sensible comes from 
the continuum limit where the Chern-Simons term (generated 
by the partons [33, 34]) leads to a mass gap [35, 36] no mat- 
ter the size of the continuum gauge coupling [37-39]. In the 
compact Abelian gauge theory, Polyakov's argument for con- 
finement [40] breaks down because the Chern-Simons term 
attaches fermion charge to the instantons; they therefore only 
contribute to fermion correlation functions, rather than to the 
free energy. 

The same physics is evident on the lattice, where the col- 
lective motion of the partons responsible for the Hall conduc- 
tivity [41] is not confined, even at large h, precisely because 
it is adiabatic. Furthermore, because the Hall conductivity 
can be expressed without reference to excitations above the 
partonic band gap [41], it is in some sense independent of 
and always contributes an order one effect (hke the continuum 
Chern-Simons term). Of course, once we move away from the 
low frequency response, we will eventually encounter the de- 
tails of the parton band structure. We note that while the lat- 
tice is essential for dealing correctly with U{1) gauge theory, 
it can effectively be dispensed with when considering non- 
Abelian groups. While the naive continuum limit of U{1) 
gauge theory fails to incorporate the instanton effects that lead 
to confinement and hence looks free, non- Abelian gauge theo- 
ries in the continuum can already accommodate confinement. 
Thus we expect that the continuum model with Chern-Simons 
term correctly captures the physics even at strong gauge cou- 
pling. 



^ We should point out that in the u = 1/2 boson case described below, the 
ring exchange term is of the same order as the interactions. 



As further evidence of the similarity between the lattice and 
the continuum, we emphasize that the deconfined phase of 
the lattice gauge theory does not have light propagating pure 
gauge excitations. Exactly as in the continuum, the strong 
gauge coupling leads to very heavy pure gauge states, namely 
the high energy "glueballs" we mentioned earlier in connec- 
tion with the physical Hilbert space. These excitations oc- 
cur at energy scale h. What we must have in the low en- 
ergy physics are renormalized partons that do not experience 
confining electric fields, as well as certain zero modes of the 
gauge field. In conjunction with the TKNN invariant, these 
gauge field zero modes lead to ground state degeneracy on 
topologically non-trivial surfaces. Finally, regarding the par- 
tons, it is useful to appeal to the bag model of confinement 
(see e.g. [42]). Imagine we start with a dilute gas of elec- 
trons where each electron may be modeled as a "bag" of the 
deconfined phase. The partons locally experience their band 
structure plus the confining potential of the bag, hence the par- 
tons propagate along the edge of the bag (a manifestation of 
the edge states in a finite size Chern insulator). When the elec- 
trons reach a density of order one, the different bags strongly 
overlap and partons may tunnel from one bag to another 
Analagous (but not identical) to the percolation transition be- 
tween quantization plateaus in the integer quantum Hall ef- 
fect, the partons are ultimately able to connect up and form a 
gapped collective state. We may picture the partons moving 
along the edges of the bags to screen out any charge sources. 
Of course, such a state does not screen in the same way as a 
metal, but a type of screening is certainly present. The sim- 
plest demonstration of this fact comes from the abelian equa- 
tions of motion, where all static fields are short ranged. 

We offer one more argument for the deconfined nature of 
the lattice gauge theory. We can consider a microscopic elec- 
tron model that possesses a completely flat Chern band at the 
cost of introducing long range hoppings. As we argue in 
Appendix A, these hoppings may be chosen to decay super- 
polynomially fast or nearly exponentially fast (although per- 
haps not exponentially fast). Now in such a model, the condi- 
tion W ^ w'^ on the electron side is vacuous, and we appear 
to be able to take h arbitrarily large on the gauge theory side. 
Yet surely such a model still realizes the fractional Chern in- 
sulator phase (at least if the partons remain gapped), since it 
closely mimics the physics of the usual continuum quantum 
fractional Hall effect. Thus we must have that the gauge the- 
ory deconfines even if h is large. We are led to conclude that 
fermions in a filled Chern band and the Chern-Simons term 
they generate have a profound impact on the dynamics of the 
gauge fields they are coupled to. 

Note that infinite range hoppings do not require a prolifer- 
ation of gauge degrees of freedom. This is because instead of 
introducing new gauge fields for each hopping, we may sim- 
ply stretch a Wilson line between the hopping sites. This Wil- 
son line only reduces the electron hopping amplitude by a fac- 
tor of the distance, a modification that can easily be absorbed 
into a redefinition of the parton hoppings without changing 
their super-polynomial decay. 
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Having dealt with the issue of confinement, we turn to the 
low energy physics. Although it is not our purpose here to 
reanalyze the low energy theory (see Refs. [3, 43]), we do 
want to make a few comments. For the Hofstadter model, 
the low energy limit of the lattice gauge theory consists 
of a non-Abelian SU{3) gauge field with a Chern-Simons 
k J [AdA + 1^3) term at level k = Cparton = 1- By level- 
rank duality, this is equivalent in the bulk to a U{1) Chern- 
Simons theory at level fc = 3 which possesses three degener- 
ate ground states on a torus. The low energy theory also con- 
tains Chern-Simons term for the background gauge electro- 
magnetic gauge field that encodes the Hall conductivity. The 
theory contains gapped fermionic matter, the partons, coupled 
to the SU{3) gauge field, and the usual anyonic statistics fol- 
low from a more detailed analysis. 

To summarize this section, we have introduced a strong 
coupling expansion that justifies our procedure for taking the 
"nth root of band structure". Furthermore, we showed that the 
unique physics of the fractional quantum Hall effect permits 
an unusually complete understanding of the connection be- 
tween the microscopic electron model and its partonic gauge 
theory description. Since the gauge theory can be reliably an- 
alyzed, our construction provides a strong argument that the 
microscopic electron model enters a fractional Chern insulat- 
ing phase. As an application of this technology, we can with 
reasonable confidence propose microscopic Hamiltonians that 
may realize non-Abelian fractional Chern insulators. We now 
turn to the case where, for definiteness, most of the attention 
has been focused - the partially filled checkerboard model - 
and describe in detail the partonic band structures and wave- 
functions generated by our approach. 

CHECKERBOARD MODEL PARTON CONSTRUCTION 

The desiderata for the nth root operation to provide a good 
fractional Chern insulator mean-field state are: 

1 . Gauge-invariant quantities should be invariant under the 
translation invariance of the original lattice model, with 
the original unit cell. This ensures that the resulting 
physical wavefunction is translation invariant. 

2. There should be a bandgap between the filled and empty 
bands. The filling fraction for each parton color is the 
same as for the original particles. 

3. The total Chern number of the filled bands should be 
nonzero. 

The final criterion (3) is crucial in order to prevent the gauge 
theory from entering a confining phase. 

In the previous section we described a strong coupling per- 
turbation theory calculation that narrows our search: define a 
new lattice with the same connectivity structure as the orig- 
inal, and define the hopping amplitude on a given link to 
be the nth root of the original amplitude. This leaves a Z„ 



phase ambiguity on each link; choosing these phases to all 
be +1 inevitably leads to partially filled parton bands; choos- 
ing these phases to break lattice translation invariance to an 
order-n subgroup produces a new bandstructure which can 
meet the above criteria. We emphasize again that the hier- 
archy of scales considered in [9-13, 31] is precisely what is 
needed to have a reliable strong coupling expansion. Thus 
our nth root procedure is, modulo the important discrete free- 
dom described above, a unique and precise mapping from the 
microscopic model to the gauge theory. For the moment we 
proceed with a general description keeping n arbitrary; later 
we will specialize to bosons (n — 2) and fermions (n — 3). 

Consider a tight-binding model defined by the following 
hamiltonian 

H = - trr'4cr'+h.C. (2) 

rr' 

Here , c denote creation and annihilation operators for elec- 
trons or hardcore bosons; the sites r are drawn from the set V 
which is a Bravais lattice plus basis. Our mean field ansatz for 
the nth root of this lattice is 

^mf = - E E {trr')'^''^rr'fUr'c.+h.C. (3) 

, fa create and annihilate n colors (labelled by a) 
of fermionic partons, from which the original electrons 
or hardcore bosons are constructed as baryons, Cr = 
^ Ea 1 . . . Q „ /ra 1 • • • /t-q „ ■ At this point, we have not yet included 
the gauge fluctuations, so the color label in the mean field 
Hamiltonian above should be treated like a flavor or band in- 
dex. We assume that the mean field Hamiltonian respects the 
full SU (n) symmetry, ujrr' are nth roots of unity expressing 
the ambiguity in the nth root operation; they can be regarded 
as variational parameters. Not all of them can be eliminated 
by rephasing the parton operators fra i^rfra i^r indepen- 
dent of a). It is important to note that although the graph is 
the same as that of the original model, the unit cell of the par- 
ton mean field Hamiltonian may be expanded relative to the 
original due to the phases ojrr' ■ 

For definiteness, we focus on the checkerboard model stud- 
ied by Sun et al [31, 44]. This is a tightbinding model with 
two sites per unit cell and hopping amplitudes out to the next- 
next-nearest-neighbor (4N). Electrons in this bandstructure at 
i filling (that is, one electron for every 3 unit cells, or every 6 
sites) with nearest-neighbor repulsive interactions, exhibit an 
insulating state [9, 10, 13] as do repulsive hardcore bosons at 
I filling [12]. The electron model has been studied for system 
sizes up to 12 electrons on 36 unit cells [13], and at various 
values of the hopping amplitudes including t" = ("7r-flux" 
model of [9]). 

These works find evidence for l/i/ degenerate groundstates 
on the torus (in each case, v is the number of particles per 
unit cell), which flow into each other under threading of flux. 
This is precisely the finite size signature of a fractional quan- 
tum Hall liquid, however, we note that an electronic charge 
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density wave which breaks lattice translation symmetry also 
shows the same qualitative physics. Indeed, in the limit of a 
thin torus, the quantum Hall fluid actually takes the form of a 
one dimensional CDW [45, 46]. However, on a torus of aspect 
ratio near one, the mixing between the three ground states are 
very different in the FCI and CDW cases; furthermore, the 
entanglement spectra will also be very different. Following 
the numerical calculations and our analytic results above, we 
assume the state is a FCI and go on to provide candidate wave- 
functions for each of these states. 

The structure of the lattice is indicated in Fig. 3 (see Fig. 1 
for a different view). We now further specialize our notation 
to the checkerboard case by explicitly distinguishing the two 
sublattices in V. More explicitly, define Fourier modes for the 
a and b sublattices (empty and filled circles in Fig. 3): 



1 



E 



1 



(4) 



where Na and N\j are the number of sites in the a and h sublat- 
tices. The Hamiltonian is H = J^keBZ i^2N + h^N + ^4w) 
with 



{kj:-\-ky) _|_ gf ( — ^a;+fct/) 



(5) 



h.c. 



J(-k^ + ky) 



Because the models we consider have fractional filling for 
the microscopic bosons and electrons, we must enlarge the 
unit cell of the partons to achieve a gap in the single par- 
ticle parton spectrum. However, we wish to insure that the 
original translation symmetry for gauge invariant variables is 
preserved (see Ref [27] for a thorough description of this re- 
quirement and its relation to the projective symmetry group). 
As an aside, we note that this is not a fundamental require- 
ment since explicit translation symmetry breaking and frac- 
tionalization can coexist, but because of its simplicity and rel- 
evance to recent numerical efforts, we choose to isolate the 
fractionalization physics. Thus we assign the discrete phases 
LOrr' in a way that doubles (bosons, n = 2) or triples (elec- 
trons, n = 3) the unit cell. We will label the subcells within 
the enlarged unit cell using script letters a, b, .... Our rule for 
assigning the phases follows: when leaving the ath 

subcell of the enlarged unit cell, particles acquire a phase w", 
where w„ = e^^*/". This has the interpretation as the sponta- 
neous development of a magnetic field of the slave gauge field, 
specifically, a field in the center of the gauge group. More ex- 
plicitly, define clock matrices in the enlarged-unit-cell basis 
a, b = l..n: 



rj(a)„b 



(6) 



We first rewrite the original Hamiltonian in the new basis, 
so far a purely cosmetic transformation. The original Hamil- 



tonian is now a sum over the reduced Brillouin zone BZ': 

H= y {a\b^)^l^,,{k){a,b)l . (7) 



E 

k£BZ' 



Now we take the nth root. Our definition of the 71th root tight- 
binding model consists of replacing f) in (7) by 

f)«throot ^ f]27V • + i)3N ' ^{a') + ^AN ' ^{a") ; (8) 

the parameters a = {a, a' , a") run from 1 to n and represent 
independent choices of nth root for each type of hopping. The 
resulting bandstructure for n — 2 and one choice of as is 
shown in Fig. 4. 




FIG. 3: The checkerboard flat band model and its cube root. The 
original unit cefl is shaded in green; the 3 x -enlarged unit cefl is less 
shaded. The arrows indicate the direction in which the hopping am- 
plitude is e"^. t'l (t'2) is the NNN amphtude associated with the solid 
(dashed) lines. 

As we have repeatedly stated, the trick is to have the 
quadratic parton Hamiltonian break the translation group up 
to gauge transformations. The simplest such Hamiltonians, 
realized in our examples below, are proportional the identity 
on the color indices. We may make use of discrete fluxes to 
obtain gapped topological bands, but these lie in the center of 
SU{n) and don't break the parton gauge group. A detailed 
analysis shows that in such cases translation invariance for 
the electrons requires the same flux (mod 2tt) through every 
equivalent loop in the enlarged unit cell. A simple way to see 
this is to note that the electron interactions, as generated by 
the strong coupling expansion, will not be translation invari- 
ant without this constraint. 

More generally, pick a parton Hamiltonian Hi that we like, 
which gives nice gapped topological bands, but whose fluxes 
do not preserve the original unit cell. Define H2 (H3) to be the 
Hamiltonian obtained by translating all the hoppings by one 
(two) original unit cell of the electrons. Now let the motion 
of the blue parton be governed by Hi, let red move accord- 
ing to H2, and let green move according to H^. The electron 
wavefunction is the product of the wave functions: ^'(z) — 
sli{z)sl2{z)sl3{z). By construction, translation by one orig- 
inal unit cell T acts by Tsli — sl2,Tsl2 — sl^jTsl^ = sli 
and therefore T"^ = ^i. 

More formally, our new hamiltonian can be described as 
follows. Our old hamiltonian was iJofd — J2a=i s^^icIcq. 
where a is the color index. The new Hamiltonian is Hnr.,,, = 
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haC'aCa- The only catch is that our new Hamiltonian no 
longer preserves the whole SU (3) parton gauge group. The 
unbroken gauge group is 

G = {U e SU{3) s.t. Hnew is preserved under — >■ U^cp}. 

Using this definition, the subgroup of SU{3) that preserves 
the new Hamiltonian is U{1) x U{1). This group acts by var- 
ious diagonal phase rotations of the parton colors. The action 
of translation by one original-unit-cell, accompanied by a per- 
mutation of the colors (the Weyl group of SU{3)) is a projec- 
tive global symmetry. The resulting Chern-Simons theory is 
[7(1)4 X U{l)e, which has the same topological properties as 
the full SUi3) theory. 

Bosons at 1/2 filling 




FIG. 4: Left: The topological flat bands found in the checker- 
board model by [31], at the optimal hopping amplitudes. Energy 
is measured in units of the NN amplitude t, and wavevectors are 
measured in units of the inverse NNN lattice spacing. The Hamil- 
tonian is the one used in the numerical work [9-13], which is mi- 
nus that studied in [31]. Right: The bandstructure for a square root 
of the checkerboard flat band model. Each parton hopping ampli- 
tude is the square root of the optimal-flatness values chosen in [31]: 

* = 1' - \/^' = ^ '4^^' ^ = 

57r/8,a = (1,1,0),^ = (0,1,1,1);. The discrete parameters S 
is explained in appendix B. The Chem number of the bottom band is 
-1. 

Since the hardcore bosons are at 1/2 filling, and the bo- 
son creation operator is 6^ = j\f\, each color of parton is 
also at 1/2 filling. Following the procedure outlined above, 
we find bandstructures like the one in Fig. 4 where the lowest 
band is separated by an energy gap from the other four The 
mean-field groundstate for each parton species is then a Slater 
determinant of the states in the lowest band. The candidate 
boson groundstate wavefunction is the projection of this state 
onto the gauge singlet sector, a procedure motivated by the 
strong coupling limit in the gauge theory. In this case, where 
the whole SU(2) gauge symmetry is preserved by the mean 
field parton Hamiltonian, the boson wavefunction is just the 
square of the Slater determinant. We see that it naturally in- 
corporates both the physics of the Chern band, since the par- 
tons have a Chern number, as well as the physics of strong 
correlation, since the boson wavefunction, as the square of a 



fermion wavefunction, forbids the bosons from approaching 
each other. 

The sum of the Chern numbers of the filled bands for each 
parton color is 1. In the continuum limit they therefore pro- 
duce a Chern-Simons term for the SU(2) gauge field with 
coefficient k = 1. The effective field theory is SU(2) level 
1 Chern-Simons gauge theory, which is related by level-rank 
duality to U{1) level 2. This state exhibits charge- 1/2 quasi- 
particles with anyonic statistics. Furthermore, each parton 
color has a Hall conductivity of (|) i, and with two col- 
ors, we have a total Hall conductivity of ^ ^ . More generally, 

an n color parton model will have a Hall conductivity of ^ ^ 
for the Abelian phases considered here. 

In the bandstructure shown in Fig. 4 we have followed very 
literally the leading-order result of the strong coupling expan- 
sion. If we permit ourselves to treat the absolute values of 
the parton hopping ampitudes as variational parameters (the 
phases are sacred), we can achieve flatter parton bands with 
larger bandgaps. This freedom should be kept in mind in fu- 
ture variational studies using these parton wavefunctions. We 
note that this bandgap matches onto the quasiparticle gap in 
the low-energy theory. 

Fermions at 1 /3 filling 

A very similar story obtains for the checkerboard flat- 
band model filled with a (spinless) electron for every six 
sites. Following the procedure detailed in the previous sec- 
tion with 71 = 3, we construct a family of cube roots of the 
checkerboard tightbinding model, parameterized by the dis- 
crete phases e^'^'"/". A favorable result is shown in Fig. 5. 
The cube of the Slater determinant of the lowest band is our 
candidate wavefunction. As before, the wavefunction incor- 
porates both strong correlation and Chern band effects. It 
is precisely analogous to Laughlin's model wavefunction for 

= 1/3 (which is also the cube of a Slater determinant). 




-1.0 -0.5 [ID 10 



FIG. 5: The bandstructure for a cube root of the checkerboard flat 
band model. Each parton hopping amplitude is the square root of 

the optimal values chosen in [31]: t — 1, t'l — ^ g+Vz ' = 
t'le-T, t" = if^^, f = 7r/12,Q = (1,1,0) . The Chern 
numbers of these bands, from bottom to top, are (1, 0, —1, 0, 0, 0). 
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CONCLUSIONS 

We have provided model wavefunctions for fractional 
Chern insulators using the parton approach. We also analyzed 
the parton gauge theory in a strong coupling expansion and 
obtained a mapping to a microscopic electron model. Further- 
more, we argued that both sides of this mapping are under 
control, based on a combination of numerical, experimental 
(in equivalent fractional quantum Hall systems), and analyti- 
cal results. Thus we have a rather special situation where the 
parton approach to fractionalized phases can reliably be used, 
even at the lattice scale. 

We find it amusing to note that the construction we have de- 
scribed is a cognate of various mechanisms for fractionaliza- 
tions of momentum [47] and D-brane charge quantum num- 
bers [48] in string theory. Like the topology of bandstructure, 
the latter enjoys a K- theory classification [49]. 

In the future, we believe it would be very interesting to 
study the gauge theory numerically, as well as to compare our 
wavefunctions to those produced in the recent numerical cal- 
culations. It would also be interesting to find a Hamiltonian 
for which our state is the exact ground state, either by push- 
ing the strong coupUng expansion or by using an analog of 
Haldane's pseudopotentials. We are currently pursuing lattice 
models that might realize non-Abelian fractional Chem insu- 
lators, and we anticipate that our strong coupling expansion 
will be useful in this endeavor. Finally, it would be quite inter- 
esting to propose analogous realistic models in 3d that might 
realize fractional topological insulators. This can be achieved 
as shown in Ref. [50] by spontaneously breaking the gauge 
symmetry in the mean field parton Hamiltonian. We can also 
consider models in 2 + 1 dimensions that break the SU(n) 
gauge group down to some subgroup, e.g. a combination of 
U(l)s or a discrete group such as Z„. These models are real- 
ized at the mean field level by parton Hamiltonians in which 
the different colors experience different hoppings, and in the 
low energy limit, they may be described by Chern-Simons- 
Higgs theories. However, as we indicated above, some care 
must be taken with lattice symmetries in such broken phases. 

More generally, our lattice gauge theory and strong cou- 
pling expansion provide powerful tools to address outstand- 
ing questions. We give three examples. First, it is experi- 
mental lore that in realistic materials the gap separating a flat 
electron band will be on the order of the spin-orbit coupling. 
Since interactions are likely to be at the Coulomb scale, it is 
potentially unlikely that the interaction energy is less than the 
electronic band gap. This brings in the danger of mixing with 
levels with different Chern number, potentially complicating 
the physics. From the perspective of our gauge theory, this 
situation can be modeled by adding bare baryon-baryon in- 
teractions. So long as these interactions are small compared 
to the bare parton gap ^ t-^, they do not appear to cause a 
phase transition, but this criterion is compatible with interac- 
tions that are larger than the bare electron gap. The story is 
only slightly complicated by the fact that the gap to quasipar- 



ticles in the gauge theory should not be set by but by some 
much smaller renormalized energy scale (one commensurate 
with electronic energy scales). Second, our mapping from the 
parton gauge theory to the electron model also suggests that 
we could have models where the bare electron band does not 
have a Chern number (despite breaking time reversal), and yet 
the partons remain in gapped Chern bands. Third, our map- 
ping naturally provides a way to take a parton band structure 
that realizes a Chern number two band and guess a micro- 
scopic electron model. This case is interesting because an 
SU{3) level 2 topological phase is non-Abelian and univer- 
sal for quantum computation. 

Since it is of particular interest, we comment a bit more 
on our progress towards non-Abelian states. Regarding the 
checkerboard model, while individual bands occasionally 
have Chern numbers with absolute value larger than 1, we 
have not found choices of a (or the other order-n ambiguities 
in the nth root procedure) that lead to bandstructures where 
the sum of the Chern numbers of the lowest bands adds up to a 
number with absolute value larger than 1. However, the model 
studied by Hatsugai and Kohmoto [7, 51] can be viewed as 
a cube root of the square lattice with NNN hoppings which 
spontaneously breaks time reversal symmetry. In a range of 
parameters, the lowest band has Chem number 2. This sug- 
gests that (spinless) electrons on this square lattice at 1/3 fill- 
ing and repulsive interactions will form a non-Abelian FQH 
state. We are in the process of investigating this claim further 
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Appendix A: Perfectly flat Chern bands 

It is possible to have a perfectly flat band with non-zero 
Chern number provided we relax the assumption of finite 
range interactions. Consider for simplicity a gapped hamil- 
tonian H with two bands with opposite and non-zero Chern 
numbers. Our construction will also work for more compli- 
cated Hamiltonians. Let P be the projector onto the lowest 
band. The hamiltonian H^^.^ = H — PHP has a perfectly 
flat lower band with the same Chem number as before. The 
question we must answer is: how non-local is the operator PI 

Consider a function f{t) whose fourier transform /(w) has 
the property that / = 1 for u in the lower band and / = 
for uj in the upper band. / is otherwise arbitrary, although we 
will want it to decay as fast as possible at large t. Since we 
require that the Fourier transform of / vanish outside a com- 
pact domain, / cannot decay exponentially fast in t (otherwise 
the Fourier transform would be analytic in a strip). However, 
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it can decay faster than any polynomial and nearly exponen- 
tially. 

Now we form the operator 



/oo 



(9) 



However, this operator is nothing but the projector P because 
of the properties of /. We are now in a position to evaluate 
{r\P\r') where rr' are positions on the lattice. The key point 
is that although we integrate over all times, long times receive 
an extremely small weighting. At early times, the amplitude 
to go from r' to r is very small because the particle hasn't had 
a chance to move. At late times, the weighting factor / is very 
small. Thus there is some intermediate time that dominates 
the matrix element, and the matrix element can be quite small 
if the separation between r and r' is large. The characteristic 
timescale in / is the inverse gap A~^. Assume for simplicity 
that the propagation under H is ballistic with speed v. Then 
the transition amplitude is very small for |r — r'| ^ vt and 
the weight factor is very small if Ai ^ 1. Combining these 
two facts, we expect a sharp decay in the matrix element of P 
once |r — r'l ^ v/A. This decay will be almost exponential 
as described above. 

Assembling everything together, the operator PHP only 
delocalizes the terms in by a size of roughly v/A up to 
corrections that are almost exponentially small. The hopping 
terms in PHP will not be of strictly finite range, but they 
will decay very rapidly. For more information about these 
techniques, see [52]. Of course, we have not proved that there 
is no finite range model with a perfectly flat Chern band. 

Regarding the gauge theory on a lattice with such long 
range hoppings, we are free to abandon our original procedure 
of associating a new SU (n) link variable with every hopping. 
Instead, we may simply stretch a Wilson line built of exist- 
ing gauge variables to connect the distant sites. This only 
modifies our strong coupling calculation by a factor of 1/x 
where x is the lattice length of the Wilson line. However, this 
modification may be reincorporated into the long range hop- 
pings without affecting their super-polynomial decay. Gauge 
kinetic terms involving these long Wilson lines will be gen- 
erated upon integrating out the fermions, but we may expect 



them to have a suitably small prefactor We can find no se- 
rious conceptual issues with the gauge theory defined with 
such long range hoppings, but we are aware that we cannot 
rule out subtle pathologies. The fact that the corresponding 
electron model appears to reliably enter a fractionalized phase 
described by a conventional low energy gauge theory suggests 
that the inclusion of these long range hoppings is not a serious 
modification of the physics. 



Appendix B: Explicit nth root Hamiltonians 

Here we give more explicit expressions for the parton hop- 
ping matrices in the checkerboard lattice models. In the 
following, the parton gauge indices are omitted because the 
Hamiltonians we describe are proportional to the identity ma- 
trix in that space. It will be convenient to define ki = kx + ky 
and fc2 = ^fc^ + ky. In each case, the tightbinding Hamilto- 
nian is a sum of nearest-neighbor (NN), next-nearest-neighbor 
(3N) and next-next-nearest-neighbor (4N) hopping terms. For 
reference, the original checkerboard model [31] is 

H = - (f)2Ar + f)3Ar + ^An) [a, bf + h.C, (10) 



BZ 



where 



e' 2 -(- e 2 



e' 2 -|- e 2 



t[ 
t' 



UN 



ifco 



1 
1 



(11) 



(12) 



(13) 



Square Root of the (upside down) Checkboard model 



To begin, define building-block matrices in the basis (ai, 61, 02, 62), where a = 1, 2 here is the subcell index: 



l2N 



( e^(T+^) 














} 



(14) 



f)3W = e}^" 



( t'2 \ 

i'l 



\ t'l / 



/ t'l \ 



t'l 

V / 



(15) 
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4N 



/O 1 0\ 

1 

10 

\ 1 / 

For n = 2, the clock matrices discussed in Eqn. (6) take the explicit form 

/ 1 \ 

10 

"laj e*"'' 

\ e*"'^ / 

Then, the general Hamiltonian for the square-root of the checkerboard model is 



(16) 



(17) 



-ff = - ^(4,61,4: 61,02,62)^ +h.c., Wnhh={h2NMNMN), n={VL{a),n{a')Ma"))- (18) 



BZ 



The parameters a, a', a" and f?i.2,3.4 take values l..n — 2 and parametrize the ambiguities in taking the square root of each 
hopping amplitude. Note that we could introduce analogs of the parameters for the 3N and 4N hoppings as well; we have not 
yet explored this possibility. 

The parameters used in Fig. 4 are 



t = i, <; = 



1 



V2' 



t'2 ~ t\e ^ , t" — it 



'^=Y' " = (1,1,0), = (0,1, 1,1) 



(19) 



Cube Root of the Checkerboard Model 

Similarly, the cube root tightbinding model is constructed as follows. 
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We introduce the VI matrices similarly: 



17(a,/3) 



/I 











^ 


1 

















e 3 

















i27rQ 

e 3 

















i2x/J 

e 3 





\0 











e 3 j 



(20) 



(21) 



(22) 



(23) 
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^ = + J2^4,blal blal 6^)(() • l1)(ai, 61, 03, ^2, 03, ^>3)^ + h.c, (24) 

with f) and f2 defined as above. This notation is slightly more general than that used in the body of the paper; the set of parameters 
we used for the cube root in Fig. 5 translates to 



t = l, i; = ^-_^, t'^=t[e^, 3^----^, <^^^, <3=(1,1,0), /3-(2,2,0). (25) 

where a = {a, a' , a"), P = (/3, /?', /3"). Note that the overall sign in (24) accounts for the sign-reversal relative to [31] which 

puts the flat band on the bottom, using the key equation (—1)'^ = — 1. 
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